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ABSTRACT 



CN ■ We present the first determination of tlie pairwise velocity dispersion (PVD) 

j^^ ' for galaxies in different luminosity intervals using the final release of the Two- 

P^ ■ Degree Field Galaxy Redshift Survey (2dFGRS). We have discovered quite sur- 

\Q • prisingly that the relative velocities of the faint galaxies at small separation are 



very high, around 700kms~^, reaching similar values as the brightest galaxies. 
O ■ At intermediate luminosities M* — 1 ( M* is the characteristic luminosity of the 

^ ■ Schechter function), the relative velocities exhibit a well defined steep minimum 

near 400kms~^. This result has been derived using a novel method to deter- 



;_i ■ mine the real space power spectrum and the PVD from the redshift space power 

^ ■ spectrum of the 2dFGRS. Combined with the observed luminosity dependence 

of clustering, our result implies that quite a fraction of faint galaxies, as well 



as the brightest ones, are in massive halos of galaxy cluster size, but most of 
the M* galaxies are in galactic halos. Our observed result is compared with the 
current halo model of galaxies of Yang et al. that was obtained by matching the 
clustering and luminosity functions of the 2dFGRS. With the model parameters 
they favored most, the halo model seems to be unable to reproduce the lumi- 
nosity dependence of the PVD because it predicts a monotonically increasing 
PVD with the luminosity. We discuss a possible solution to this model by raising 
the faint end slope of the conditional luminosity function in rich clusters. The 
PVD luminosity dependence may also be an important constraint in general on 
theories of galaxy formation, such as semi-analytical models and hydro/N-body 
simulations of galaxy formation. 

Subject headings: galaxies: clustering - galaxies: distances and redshifts - large- 
scale structure of Universe - cosmology: theory - dark matter 
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1. Introduction 

The clustering of galaxies in the Universe is characterized by their spatial positions, and 
by their peculiar velocities that lead to deviations of their motion from the pure Hubble flow. 
The big redshift surveys assembled in recent years by the diligent work of many astronomers 
give angular positions and redshifts for large numbers of galaxies. A rough 3D map can be 
obtained by placing the galaxies at distances along the line of sight derived via Hubble's law 
from their redshifts. The peculiar velocity, however, also contributes to the redshift, and this 
leads to a misplacement of the galaxy away from its true location. The local gravitational 
field is the cause of the peculiar motion, thus the redshift distortion in the galaxy maps can 
give information on the underlying matter distribution. 

The amplitude of the distortions can be estimated from the pair distribution of galax- 
ies. For pairs of galaxies at distances much larger than their separation, one can use a 
plane-parallel approximation. In the linear approximation (or large separation), the power 
spectrum of their distribution in redshift space is (Kaiser 1987), 

P^(k) = (l + /3/x2)2p(A;). (1) 

Here n is the cosine of the angle between the wave vector and the line of sight. The linear 
redshift distortion parameter (3 that is related to the linear growth factor f{^o) — f2°'^ (f2o 
is the matter density) and the linear bias factor b of the galaxies hy (3 = Q^^/b , can be 
estimated if P'^(k) can be measured on a sufficiently large scale. P{k) is the power spectrum 
in real space. The dependence on b expresses the fact that there is a bias in the galaxy 
distribution, i.e. there are differences between the galaxy and the dark matter distributions. 
The linear bias relation is just a constant ratio of the galaxy C, and the dark matter C,dm 
two-point correlation functions (2PCF) ^ = b'^^dm- It has been shown that Eq.(l) is valid 
only on sufficiently large scales (perhaps ^20/i~^Mpc or k < 0.1/iMpc"^, Scoccimarro 
2004). On smaller scales, the virial motion of galaxies within groups and clusters contributes 
significantly to the deviation of the redshift distortion from the prediction of Eq.(l). In fact, 
on small scales (< 5 /i~^Mpc), the virial motion of galaxies dominates the redshift distortion. 

Thus the redshift distortions on small scales can be used to determine the pairwise 
peculiar velocity dispersion (PVD) of galaxies. Assuming certain functional forms for the 
distribution function of the pairwise velocity (say, an exponential form, Peebles 1976) and 
for the average infall velocity, a model can be constructed for the redshift 2PCF which 
approximates the real situation well, when the coupling between the peculiar velocity and 
the spatial density of the galaxies is weak (Fisher et al. 1994b; Juszkiewicz, Fisher, &Szapudi 
1998; Davis & Peebles 1983, hereafter DP83). A comparison of the model with observations of 
the redshift 2PCF provides a test for the validity of the assumptions (such as the distribution 
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function of the pairwise velocity) and a determination of the PVD 0-12 (r). DP83 apphed 
this method to the CfA redshift survey, and determined the PVD to be 340 ± 40kms~^ at 
projected separation r^ = 1 h~^Mpc. Based on an extensive study of the PVD for all redshift 
surveys (typically containing 2000 galaxies) available before 1993, Mo, Jing, & Borner (1993) 
showed that the PVD measured is very sensitive to the presence of rich clusters in the survey. 
Because the CfA survey used by DP83 is too small to fairly represent the population of rich 
clusters in the Universe, Mo et al. pointed out that the PVD estimate given by DP83 is 
likely to be significantly biased low based on their analysis of the PVD for different surveys. 
Subsequent analyses for CfA2 by Marzke et al. (1995) and for CfA by Somerville, Davis, & 
Primack (1997) have confirmed the conclusion of Mo et al.. 

With the Las Campanas redshift survey (LCRS) of galaxies, Jing, Mo, & Borner (1998, 
hereafter JMB98) made the first accurate determination of the PVD using the above method. 
They used mock catalogs generated from N-body simulations of the concordance A cold dark 
matter (LCDM) model to correct for the observational effects including the fiber-collisions, 
and to make a fair estimation of errors for the measured PVD. They demonstrated that the 
errors of their measured PVD are between 50 and 100 kms^^, and claimed that the PVD can 
be accurately measured with the LCRS, in contrast with all previous work based on smaller 
surveys. The PVD at the projected separation 1 h~^Mpc is 570 ± 80kms~^. This result has 
been verified by Zehavi et al. (2002) with the early data release of the Sloan Digital Sky 
Survey (SDSS). For IRAS galaxies, the PVD is lower than that of optical galaxies (Fisher 
et al. 1994a,b; Mo, Jing, & Borner 1993), and the PVD is very low on small scales (Jing, 
Borner, & Suto 2002; Hamilton & Tegmark 2002). These results of IRAS galaxies are also 
consistent with the SDSS study (Zehavi et al. 2002) where it was found that the PVD of 
the blue (young) galaxies is as low as that of IRAS galaxies. Cosmological hydrodynamic 
simulations of galaxy formation seem to lead to similar conclusions with regard to the PVD 
of different galaxy populations (Pearce et al. 2001; Yoshikawa, Jing, & Borner 2003; Berlind 
et al. 2003; Weinberg et al. 2004). 

Another way to measure the PVD is to use the redshift space power spectrum. The 
redshift power spectrum P^{k,fi) can be written as (Peacock & Dodds 1994; Cole, Fisher, 
& Weinberg 1995) 

P^ik.fi) = P{k){l + ^f?fD{k^a^2{k)) (2) 

where the first term is the Kaiser linear compression, and the second term D is the damping 
effect caused by the random motion of the galaxies. For the exponential distribution function 
of the pairwise velocity, the function D is Lorentz form 
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Jing & Borner (2001a) have studied the redshift power spectrum in typical CDM models for 
different tracers -the primordial density peaks model, the cluster-underweighted scheme of 
particles, and pure dark matter, and found that different tracers may have different forms 
of the damping function. Thus, the damping function can serve as a constraint on galaxy 
formation models. From the analysis of the LCRS catalog, both Landy, Szalay, & Broadhurst 
(1998) and Jing & Borner (2001b) found that the damping function of observed galaxies is 
very close to Eq.(3). Adopting Eq.(3) for the damping function, Jing & Borner (2001b) 
measured the PVD cruik) for LCRS by setting f3 = 0.45 and found that the measured 
PVD is consistent with the PVD that JMB98 reported on the redshift 2PCF measurement. 
Hawkins et al. (2003) used Eq.(2) and Eq.(3) to model the redshift 2PCF with the Fourier 
transform, and measured the /3 and o"i2 (assuming cruik) = const.) for the Two-degree Field 
Galaxy Redshift Survey (2dFGRS) by best fitting the observed 2PCF. The advantage of 
using the redshift power spectrum to determine the PVD is that it is simple and accurate 
to model the infall effect. 

Here we use the data of the final release of the 2dFGRS^ for such an analysis to study 
the luminosity dependence of the PVD. The parameter f3 could, however, not be determined 
independently for each luminosity subsample, since there are still large statistical fluctuations 
on large scales fca/i < 1. Therefore, we fix /5 at a reasonable value of 0.45. Later we come 
back and consider the influence of a luminosity dependence (as in Norberg et al. 2002a) of 
(3. Inspection of Fig. (9) shows that the results are basically unchanged, implying that the 
results are robust to reasonable changes of the jS values. 

The 2dFGRS has already been analysed statistically with respect to the PVD (Hawkins 
et al. 2003) and the determination oi f3 (Peacock et al. 2001; Hawkins et al. 2003). But here 
we use a novel method to estimate the PVD. Furthermore, we take advantage of the large 
number of galaxies in the 2dFGRS to bin the galaxies in different luminosity intervals, and 
study the luminosity dependence of the PVD. This was not possible up to now, and we shall 
see that remarkable tests of galaxy formation models become possible with the statistical 
results presented here. The luminosity dependence of the clustering (the 2PCF) of galaxies 
in the 2dFGRS has been investigated (Norberg et al. 2001, 2002a,b) but not of the PVD. 

It has been demonstrated already in an analysis of the LCRS ( JMB98) that the observa- 
tionally measured 2PCF and the PVD require a scale dependent bias model in order to find 
acceptable theoretical fits. The so-called cluster-underweighted model that was invented by 
JMB98 to bring a flat Q = 0.2 CDM model into agreement with the data has meanwhile 
been developed into a more sophisticated "halo-occupation-distribution" (HOD) model (e.g. 



^Available at http://www.mso.anu.cdu.au/2dFGRS 
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Sheth et al. 2001; Berlind & Weinberg 2002; Kang et al. 2002). In particular, Yang, Mo, & 
van den Bosch (2003, hereafter Y03) introduced the conditional luminosity function $(L|M) 
to characterize the luminosity distribution of galaxies within a halo of mass M, and mea- 
sured the parametrized form of $(L|M) by matching the observed luminosity function of 
galaxies and the luminosity dependence of galaxy clustering in the 2dFGRS. This model has 
been applied to predict the clustering and velocity statistics of galaxies in the DEEP2 survey 
(Yan, Madgwick, & White 2003). By construction, the model of Y03 is able to reproduce the 
luminosity function and the two-point correlation function of galaxies. But the PVD, espe- 
cially its luminosity dependence, can serve as an independent test on the model, as the PVD 
on small scales is rather sensitive to the HOD. We will show that the model most favored by 
Y03 is unable to match our measured luminosity dependence of the PVD, indicating that 
some of the parameterizations adopted by Y03 have to be amended. A comparison with the 
semi-analytic models for galaxy formation (Diaferio et al. 1999; Benson et al. 2000) will be 
presented in a subsequent investigation. 

In the following chapters we describe our way of using the redshift power spectrum, the 
construction of mock catalogs from our simulations, tests of the method, the results, and the 
comparison with the HOD model of Y03. In the final discussion we summarize our results, 
and suggest construction principles for better models. 



2. Observational sample and random sample 

We select data for our analysis from the final release of the 2dFGRS (Colless et al. 2003, 
2001, hereafter COl). The survey covers two declination strips, one in the Southern Galactic 
Pole (SGP) and another in the Northern Galactic Pole(NGP), and 99 random fields in the 
southern galactic cap. Each redshift determination is assigned a quality class Q in the 5- 
point system according to the measurement accuracy based on emission and absorption lines. 
Quality Q = 1 or 2 means a doubtful redshift estimate, (5 = 3 means a probable redshift 
with the confidence 90 percent, and Q = 4 or 5 means a reliable redshift. Quality classes 1 
and 2 are considered failures. The redshift sampling completeness R{0) {9 is a sky position) 
that is defined as the fraction of targeted galaxies for which a redshift is measured with 
Q > 3, is available for each sky sector (COl). In this paper, only galaxies in the two strips 
are considered. Further criteria for galaxies to be included in our analysis are that they are 
within the redshift range of 0.02 < z < 0.25, have the redshift measurement quality Q > 3, 
and are in the regions with the redshift sampling completeness R{9) better than 0.1 (where 9 
is a sky position). The redshift range restriction should ensure that the clustering statistics 
are less affected by the galaxies in the local supercluster, and by the sparse sampling at 
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high redshift. The redshift quality restriction is imposed so that only galaxies with reliable 
redshift are used in our analysis. An additional reason is that the redshift completeness 
mask provided by the survey team, which is used in our analysis, is constructed for the 
redshift catalog oi Q > 3. The last restriction is imposed in order to eliminate galaxies in 
the fields for which the field redshift completeness cp is less than 70 percent (see COl about 
the difference between R{0) and cp)- These fields have not been included in computing the 
redshift mask map R{6). As a result, there are a total of 190504 galaxies satisfying our 
selection criteria, 78190 in the NGP strip and 112314 in the SGP strip. 

In order to study the luminosity dependence of the PVD, we divide the galaxies into 
10 subsamples according to their absolute magnitude. The subsamples are successively 
brightened by 0.5 magnitude from the faintest sample Mj, = —17.0 + 51og]^Q/i to Mf, = 
—21.5 + 51og]^Q/i, with successive subsamples overlapping by 0.5 magnitude. Here h is the 
Hubble constant in units of lOOkms^^Mpc^^. The details of the subsamples studied in this 
paper are given in Table 1. For computing the absolute magnitude, we have used the k- 
correction and luminosity evolution model of Norberg et al. (2002b, k + e model), i.e., the 
absolute magnitude is in the rest frame bj band at z = 0. We assume a cosmo logical model 
with the density parameter Qq = 0.3 and the cosmological constant Aq = 0.7 throughout 
this paper. 

A detailed account of the observational selection effects is released with the catalog by 
the survey team (COl). The limiting magnitude changes slightly across the survey region 
due to further magnitude calibrations that were carried out after the target galaxies had 
been selected for the redshift measurement. This observational effect is documented in the 
magnitude limit mask b^j^{6) (COl). The redshift sampling is far from uniform within the 
survey region, and this selection effect is given by the redshift completeness mask R{6). The 
redshift measurement success rate also depends on the brightness of galaxies, making fainter 
galaxies more incomplete in the redshift measurement. The fi{6) mask provided by the survey 
team is aimed to account for the brightness-dependent incompleteness. These effects can 
be corrected in computing the redshift two-point correlation function through constructing 
random samples that properly include these selection effects. We generate random samples 
in the same way as described in Jing & Borner (2004). Each random sample for a northern 
or southern luminosity subsample contains 100,000 random points. 



3. Redshift power spectrum and the pairwise velocity dispersion 

We measure the redshift two-point correlation functions ^2(s) following the method of 
JMB98. 
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We convert the redshift two-point correlation function to the redshift power spectrum 
by the Fourier transformation: 

P^(k) = y"e.(s)e-^''-^rfs. (4) 

In cyhndrical polar coordinates {vp, 0, tt) with the 7r-axis parallel to the line-of-sight, P'^(k) 
depends on kp, the wavenumber perpendicular to the line-of-sight, and on kj,-, the wavenumber 
parallel to the line-of-sight. The power spectrum can be written 

ps^kp, k^) = I i,{rp, 7r)e-^[^^"^^°^(^)+"^-Vprfrprf0rf7r . (5) 

With some elementary mathematical manipulation, we get the following expression: 

/»oo /»oo 

P^ {kp, k^) = 27r d-K I rpdrp^^{rp,TT)cos{k^7r)Jo{kprp) (6) 



where Jo{kprp) is the zeroth-order Bessel function (Jing & Borner 2001b). 

^z{i^p, tt) is measured in equal logarithmic bins of Vp and in equal linear bins of n. The 
reason why different types of bins are chosen for Vp and tt is the fact that C,z{i^p, tt) decreases 
rapidly with Vp but is flat with vr on small scales. Thus this way of presenting C,z{i^p,'^) is 
better than using the log-log or the linear-linear bins for Vp and tt, and is also suitable for the 
present work. The peculiar velocity of a few hundred kms~^ should smoothen out structures 
on a few h~^Mpc in the radial direction, and the linear bin of Attj = 1 h^^Mpc is suitable 
for resolving the structures in the radial direction. With logarithmic bins chosen for r^, the 
Tp dependence is resolved well, because otherwise the small scale clustering on the projected 
direction cannot be recovered. With this bin method, we obtain the power spectrum: 

P^{kp, k^) = 27r ^ AvTirp jA lnrpj^^{rpj, ttj) cos{k^7ri)Jo{kprpj) (7) 

where ttj runs from —50 to 50 h^^Mpc with Attj = 1 h^^Mpc and r^j from 0.1 to 50 h^^Mpc 
with Alnrpj = 0.23 (Be careful not to confuse two tts in Eqs. (6) and (7): the flrst n in the 
right-hand-side has the conventional meaning, i.e. 3.14159..., and the others are for the axis 
along the line-of-sight.). We make the summation of Eq.(7) with rectangular boundaries in 
vr and Vp. 

The fluctuations of ^zi'f^p,'^) at large separations may bring errors to the determination 
of the redshift power spectrum. We improve the measurement by down weighting C,z{i^p,tt) 
at the larger scales. Speciflcally, we use a Gaussian window function 

W,is) = expi-^^) (8) 



to weight the two-point correlation function with S = 20h ^Mpc. That is, the measured 
redshift power spectrum is 

P^'^^ikp, K) = 27r ^ ATVivljA \nrpj^^{rpj, Tri)Wg{rpj, ttj) cos{K7Vi)Jo{kprpj) . (9) 

id 

The weighting reduces the noise in the measurement, but introduces a systematic bias 
to the redshift power spectrum. This is the reason why we distinguish this measured redshift 
power spectrum P^'"^{kp, k^^) from the expression P'^{kp, k^^) in Eq.(7). But they are related 
by the following equation, 

^''•'"(k) = -^ J P'ik,)Wgik-k^;S)dk^ (10) 



and 



W,i^-'S) = j^^±^^eM-ls'k^)- (11) 



The weighting may change the redshift power spectrum at A; < n/S = 0.15 /iMpc"^, 
making the spectrum more isotropic and biasing the value of P''^{kp, kj^) generally. This effect 
can be estimated as follows. We take Eqs.(2) and (3) for the redshift power spectrum P'^(k). 
As a typical example, we take jS = 0.45, a^ = 500kms~^, and a linear Cold Dark Matter 
power spectrum from Bardeen et al. (1986) with the shape parameter F = 0.2 for P'^(k), 
and compute P'^'"^(k) with Eq.(lO). The result is plotted in Fig. (1) which is also compared 
with P'^(k). The weighting makes the redshift power spectrum significantly rounder at A; = 
0.1/iMpc~^, but the effect becomes negligible for k > 0.2/iMpc^^. Furthermore, the figure 
shows that the P^(k, fi = 0) at /i = is changed little even at A; = 0.1 ZiMpc"^, indicating 
that the real space power spectrum may be measured unbiasedly for A; > 0.1 hMpc~^ (because 
P(A;) = p3{k,fi = 0)). 

A further test of the method is carried out by comparing the redshift power spectrum 
of a simulated galaxy catalog. The galaxy catalog is generated using the halo model (Y03, 
Yang et al. 2004, hereafter Y04) for a LCDM simulation of boxsize 300 h'^Mpc. The details 
about the CDM model, the simulation, and the method of generating the galaxies will be 
given in the next section. For our test, we select galaxies with luminosity in —19.5 < 
Mf, — log^o h < —18.5. We computed the redshift two-point correlation function by counting 
the galaxy pairs, and measure the redshift power spectrum using the procedure outlined 
above. These spectra are compared with the redshift power spectra measured with the 
Fourier transformation (FT) of the galaxy density field. This is shown in Fig. (2), with the 
symbols for P{k,fi) obtained with FT, and the line for those measured through the 2PCF. 
The power spectra measured with FT are free of the weighting effect [Eq.(lO)]. The figure 



-9- 



confirms the results of Fig. (1), and indicates that our method can give an unbiased estimate 
of P{k,ji) for k > 0.2/iMpc~^ from which we will measure the velocity dispersion (T^{k) 
through equation (2). The real space power spectrum P{k), which is mainly determined 
from P^{k, /i = 0), can be determined for A; > 0.1 /iMpc~^. 



4. Simulation, Halo model, and mock catalogs 

We simulate galaxy catalogs using our cosmological N-body simulations. The cosmologi- 
cal model considered is a currently popular flat low-density model with the density parameter 
Qq = 0.3 and the cosmological constant Aq = 0.7 (LCDM). The shape parameter F = flQh 
and the amplitude ag of the linear density power spectrum are 0.2 and 0.9 respectively. 
Two sets of simulations, with boxsizes L = 100 /i~^Mpc and L = 300 /i~^Mpc, that were 
generated with our vectorized-parallel P^M code (Jing & Suto 2002; Jing 2002), are used in 
this paper. Both simulations use 512^ particles, so the particle mass rrip is 6.2 x IO^H^^Mq 
and 1.7 x IO^^H^^Mq respectively in these two cases. We have four independent realizations 
for each boxsize. Dark matter halos are identified with the friends-of-friends method (FOF) 
using a linking length b equal to 0.2 times of the mean particle separation. All halos with 
ten members or more are included for generating the galaxy catalog. 

We populate the halos with galaxies according to the prescription proposed by Y03. 
The luminosity function of galaxies in a halo is assumed to be a function of the halo mass 
M, $(L|M) that is further parametrized as the Schechter function. The parameters of the 
Schechter function are functions of M, further parametrized through the mass-to-light ratios 
of halos. There are a total of 8 parameters (without classifying galaxies with spectral types) 
that are determined by best fitting the observed luminosity function and the luminosity 
dependent clustering of galaxies in the 2dFGRS. It is shown by Yang et al. that these two 
observations are well reproduced by their model, but there is some degeneracy in the model 
parameters. In our work, we adopt model Ml in Y03 for the parameters of the halo model, 
and populate the halos with galaxies (the luminosity, position and velocity) in a similar 
way as Y04. The code for populating the halos with galaxies was, however, written by 
us independently. We adopt the "FOF satellites" scheme of Y04. We assign the satellite 
galaxies of a halo the position and velocity of dark matter particles randomly selected from 
the halo. Instead of locating the central galaxy at the center-of-mass of a halo as in Y04, 
we locate it at the potential minimum of the halo and assign its velocity with the halo 
bulk velocity. The mass resolution of the simulations sets a faint limit for which galaxies 
can be presented. According to Y04, galaxies are complete in the 300 h~^Mpc simulation 
for magnitude brighter than Mb — Slog^g^ — —18-0 and in the 100/i~^Mpc simulation for 
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magnitude brighter than Mh — Slogj^g/i = —15.0. In Fig. (3), we present the luminosity 
functions of the model galaxies which do agree well with the observation of the 2dFGRS, 
and also agree with each other between the two sets of simulations. 

It is straightforward to produce mock catalogs of the 2dFGRS, since we have the infor- 
mation of positions, velocity and luminosity for each galaxy. We first duplicate the simulation 
volumes periodically along the main axes, and execute the selection effects according to the 
observational masks. We produce 5 mock galaxy catalogs for each realization of a simu- 
lation, so we have a total 40 mock samples. For studying the redshift power spectrum of 
galaxies fainter than Mfe — 5 log^Q h = —18.0, we will use the mock samples of the 100 h~^Mpc 
simulations, otherwise we will use those of the 300 h~^Mpc simulations. 

It is important to understand the numerical artifacts of the simulations that may affect 
the result here. We believe that the force resolution of the simulations that adopt P^M, are 
sufficient for the current study, since only the global information of halos is needed. The 
mass resolution has already been discussed above. It can be easily taken into account in our 
analysis. The boxsize may be the only concern that we need to consider seriously, because 
the most massive halos may be under-represented for the lack of long-wavelength density 
fluctuations, and the redshift power spectrum cannot be measured around the fundamental 
wavelength of the simulation. Since we measure P{k, fi) for A; > 0.1 /iMpc^^, we believe that 
the boxsize 300 h~^Mpc is sufficiently large both for studying the large-scale clustering and 
for fairly sampling the massive clusters. To quantify the possible effect of the 100 h~^Mpc 
boxsize, we measure the redshift power spectrum for the luminosity —19.0 < M^- 5 log^g ^ < 
— 18.5 from the mock catalogs. The results of the two sets of mock catalogs are compared 
in Fig. (4), which clearly show that the simulation of 100/i~^Mpc box size is sufficiently 
big in volume for studying the redshift power spectrum for k > 0.16 /iMpc^^. Also the 
population of the halos that host the faint galaxies in the simulation of 100 h^^Mpc should 
not be affected by the limited volume. 



5. Results 

The statistical quantities are defined in the Fourier space (k) in the present work, which 
are often compared those obtained in the coordinate space (r) in previous work. When 
doing this comparison, one uncertainty is the scale correspondence between the k and r 
spaces. Although the correspondence is r = 27r/A; in mathematics, we think that it is more 
appropriate to compare the two-point clustering and the PVD with the relation r <^ 1/k 
based on the following facts. The two-point clustering can be expressed as a sum of two- 
halo and one-halo contributions (e.g. Kang et al. 2002). The transition from the two-halo 
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term dominance at large scales to the one-halo term dominance happens at ~ 0.5 /iMpc~^ 
in /c-space or at r ~ 2 h~^Mpc in the r-space for the concordance LCDM model(e.g. Kang et 
al. 2002, Y04). When comparing the scale dependence of PVD, Jing & Borner (2001a) also 
argued for k <^ 1/r. Therefore, we adopt r <^ 1/k in the current paper when comparing the 
statistical quantities in the two spaces. 



5.1. The luminosity dependence of the PVD 

Our main results are shown in Figs. (5) to (8). In the introduction we have said that 
it should be possible to derive f3, 0"i2, and P{k) from a measurement of the redshift space 
power spectrum P'^{k,fi). However, because we measure P^{k,fi) only up to the scale 
k = 0.1 hMpc~^, there exists a strong degeneracy in determining the parameters 0"i2 and P 
(Peacock et al. 2001) from P^{k,fi). Moreover, au could be a function of A;(Jing & Borner 
2001a). We therefore fix /9 = 0.45 as a reasonable estimate (Hawkins et al. 2003), and then 
determine P{k) and a^ from the data. There could also be some luminosity dependence 
in (3, and we investigate this by using the luminosity dependence of the bias parameter 
b/b* = 0.8 + 0.2(L/L*) given in Norberg et al. (2002a) and j3* = 0.45, where the quantities 
with superscript * are those at the characteristic luminosity M*. The result is shown in 
Fig. (9). The PVD for the brightest galaxies is slightly increased compared to the result in 
Fig, (8), and it is slightly decreased for the faint galaxies. The changes are small, therefore 
it seems justified to use a constant p. In Figs. (8) and (9) we have also plotted the PVD 
values at the number-weighted luminosity value in each bin. 

In Fig. (5) the basic measurement of the power spectrum in redshift space P^{k,fi) 
is shown. The four panels in this figure correspond to four different luminosity intervals 
(samples 3,5,7, and 9 of table 1) from faint to bright galaxies. The symbols are the results 
for the full 2dFGRS, the dotted lines for the south, and the dashed lines for the north sample 
for each A;-value. The values of k range from 0.20 /iMpc^^ at the top to 3.2 /iMpc^^ at the 
bottom with an increment of A log^^Q k = 0.2. The south and north samples agree quite well 
with the full survey indicating that cosmic variance is not a problem. The solid lines are 
the best fits obtained by applying equation (2) to the data of the whole survey. The power 
spectrum for the larger A;-values decreases quite strongly with n, more than a factor of 10 
between /i = and fi = 1. There is a small luminosity dependence. P'^{k,fi) increases for 
the bright samples by a factor ~ 2 for all A;. A similar dependence on luminosity for the 
2PCF was found by Norberg et al. (2002a) (see also the real space spectrum below). 

In Fig. (6) the real space power spectrum of the 2dFGRS P{k) is displayed for the four 
different luminosity intervals (symbols). Again we can see that P{k) for brighter samples 
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are higher. The error bars are derived from the mock samples described in the previous 
section. We prefer this error estimate to the bootstrap method which in general gives 50% 
smaller errors. The mock sample errors are just the standard deviations occurring, when 
the identical analysis is carried out for the 20 mock samples. These error bars adequately 
include the sample to sample variations (the cosmic variance), and they are not sensitive 
to the bias model used (JMB98). The figure shows that P{k) is quite reliably determined 
for k between O.l/iMpc"^ and 4/iMpc~^. Although P{k) = P^{k,jj, = 0) in principle, 
P{k,fi = 0) fluctuates around the true value of P{k) in practice, because in a finite sized 
survey the number of Fourier modes is always limited. P{k) is better determined, if P^{k, fj,) 
at different angles is combined with P{k) being treated as a free parameter, as is done here. 
The spectrum is approximately a power law for the range of k considered here. It decays with 
k approximately as k~^-^, with the slope of the brightest sample somewhat shallower. The 
smooth solid lines are the predictions of the halo model of Yang et al. which is implemented 
as described in section 4. The agreement with the data is satisfactory on large scales where 
k < 0.5 /iMpc"^, not a surprising fact, since the halo model has been constructed such as to 
reproduce the luminosity-dependent clustering length of the 2dFGRS (i.e. linear to quasi- 
linear regimes). Nevertheless, we may note that for samples 5,7, and 9 the halo model gives 
an indication of a change in slope of P{k) at k = 1/iMpc^^. This can be understood as 
the transition between the scales where the pair counts are dominated by galaxy pairs in 
the same halo to the larger scales, where pairs of galaxies mostly are in separate halos. In 
the data such a change in slope is much less pronounced, but it is present for the brightest 
samples (see also Y04, Zehavi et al. 2003). The discrepancy between the model and the 
2dFGRS data at small scales of A; ^ 1 /iMpc~^ is very significant, at ~ 3a level for sample 9, 
at ~ lOo" level for sample 7, and at ~ 5o" level for sample 5. We also note that the clustering 
of the faint sample (sample 3) is higher than the halo model on all scales at ~ 1.5o" level. 

The PVD is measured simultaneously with P{k), and the results are presented in 
Fig. (7). It seems that for the /c- values used here, (Ji2{k) is a well-determined quantity. 
For the faintest sample, the PVD is persistently higher than that of the halo model (at 
~ 2o" level), indicating that the faint galaxies in the Universe may reside in more massive 
halos than the halo model predicts. This can also explain why the P{k) of this sample is 
systematically higher than the halo model prediction. For sample 5, the difference in the 
PVD between the 2dFGRS and the halo model is only 50 to lOOkms"^, consistent with the 
fact that their P{k) are also close, though the agreement is not very good as the discrepancy 
between the model and the observation is significant at 2a level. For the bright galaxies 
(sample 7 and sample 9), the halo model prediction rises quickly with k at k > l/iMpc~^, 
as the galaxy pairs on these scales are mainly contributed by those from the same massive 
halos. This behavior was also seen in the power spectrum P{k) at small scale where P{k) 
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rises quickly with k. Especially for sample 7, the observed PVD does not rise with k and 
is only around 500kms~^, which is expected if the galaxies are more distributed among 
isolated halos with the velocity dispersion of 200 ~ 300kms~^ (bright spirals). For sample 
9, the PVD rises with k, but much slower than the halo model around 1 hMpc"^, which also 
indicates that the brightest galaxies in the real Universe are distributed more in halos of 
different mass than the halo model predicts. The difference between the halo model and the 
observation is significant at > 3a level for the bright samples. 

The luminosity dependence of the PVD is shown most clearly in Fig. (8), where we 
have plotted ai2 at k = 1 /iMpc~^ for the 10 overlapping samples listed in table 1. The 
surprising result is a strong dependence on luminosity with the bright and the faint galaxies 
reaching high values of ~ 700kms~^, and with a well defined minimum of ^ 400kms~^ for 
the galaxies of sample 8. The bright and the faint galaxies apparently have high random 
motions, as expected for objects in a massive halo or in a cluster. The M* like galaxies are 
rather moderate in their PVD, and probably reside in galaxy size halos. The thick solid line 
represents the prediction based on the halo model of Y04 (§4) that clearly does not match 
our observation of the PVD. 



5.2. Implications for the halo model 

The results of Y04 as well as of the current work show that the halo model adopted in 
this paper can reproduce the luminosity function very well, but cannot match the clustering 
data of galaxies perfectly. The model can quite successfully account for the luminosity 
dependence of clustering on large scales (i.e. the luminosity dependence of the clustering 
length), though it seems to underpredict for the faintest galaxies (Figure 6). Moreover, the 
clustering at small scales (r < 2 h~^Mpc or k > 0.5 /iMpc^^) predicted by the model is 
higher than the 2dFGRS data by a factor of 2 to 3. Also the PVD and the quadrupole- 
to-monopole ratio of the redshift space correlation function are higher in the model than 
the 2dFGRS data when the luminosity classification has not been applied. To solve these 
problems (or part of them), Y04 proposed several possible solutions, that is, to raise the 
mass-to-light ratio to 1OOO(M/L)0 for rich clusters in the B-band, to introduce a strong 
velocity bias 6vci = o'gai/o'dm = 0.6 (where Ugai and /(Tdm are the velocity dispersions of 
galaxies and dark matter in halos), or to lower the clustering amplitude as to about 0.7. 
As Y04 realized, the first two solutions are not realistic, because the observed mass-to-light 
ratio is (450 ± 100) (M/L)© (Fukugita et al. 1998) and (363 ± 65) (M/L)© (Carlberg et al. 
1997) for rich clusters in the B-band, and hydro/N-body simulations show that velocity bias 
is rather minor for rich galaxy systems (Yoshikawa, Jing, & Borner 2003; Berlind et al. 2003). 
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Therefore we only consider the last possibility to see if a low as = 0.7 can bring the halo 
model into agreement with our luminosity-dependent PVD result. 

In our N-body data library, we have one simulation with L = 100 h^^Mpc and another 
with L = 300 h~^Mpc in a ag = 0.7 LCDM model. The model and simulation parameters of 
these two simulations are the same as those of the simulations with the same boxsize L in §4, 
except for that the density parameter Qq = 0.25 is slightly lower for the new model. Xiaohu 
Yang has kindly provided us with the parameters of the halo model for this cosmological 
model. We find that the luminosity function of the 2dFGRS can be reproduced and the 
clustering is improved on small scales (but still higher for luminosity bins around M*) in 
this model. Its PVD is shown in Fig. 8 as a function of luminosity, which is considerably 
lower than the data of 2dFGRS. Compared with the low ag, model of Y04, our Qq is slightly 
lower. Since the PVD is approximately proportional to as^Q^, we can boost the PVD by 
(0.30/0.25)*''® which is also shown in Fig.8. We can see that this PVD after correcting for 
the Qq difference is still lower than our 2dFGRS data. More interestingly, this low ctq model 
also predicts a monotonically decreasing PVD with the decrease of the luminosity, implying 
that simply changing the value of as cannot overcome the difficulty of the halo model to 
explain the bimodal nature of the luminosity dependence of the PVD. 

It should be possible to use the PVD to constrain the value of as if the shape of the 
PVD is reproduced in the halo model. Y04 argued for as = 0.7 from a comparison with the 
2dFGRS PVD without a luminosity classification (Hawkins et al. 2003) (see also van den 
Bosch, Mo, & Yang 2003, for other independent arguments), but our comparison of the halo 
model with our 2dFGRS data seems to favor a higher as, since the PVD of the this low 
as model is lower than the observed data at all luminosities. Considering that most of the 
2dFGRS galaxies are in the luminosity interval —20.5 < Mf, — 51og]^Q h < —18.5, we would 
require as ~ 0.82 for the Qq = 0.3 according to the as^^'^ scaling if we want the halo model 
to match the observed PVD at Mf, — 5 log^g ^ — ^19.5. This value of as is in good agreement 
with the recent determinations based on the WMAP and the SDSS galaxy-galaxy lensing 
data (Spergel et al. 2003; Seljak et al. 2004a). Although this is not a rigorous determination 
for as because we have not yet reproduced the shape of the luminosity dependence of PVD, 
our comparison indicates that the concordance model favored by the WMAP and SDSS data 
is consistent with our PVD data. 

The PVD is an indicator of the depth of the local gravitational potential. Therefore we 
are inevitably led to the interesting conclusion that the bright and the faint galaxies move 
in the strongest gravitational field. A substantial fraction of them must be in clusters, while 
most galaxies around magnitude M^ — blogi^h = —20.5, the M* galaxies in the Schechter 
luminosity function, populate the field. We believe this result, that is, the shape of the PVD 
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luminosity dependence (not its absolute amplitude), constitutes a new challenge to the halo 
model of Yang et al.. The straightforward implication for the halo model is to increase the 
faint galaxy population in rich clusters. With the parameters of the Yang et al.'s halo model, 
the faint end slope 0:15 of their conditional luminosity function in rich clusters is one of the 
crucial quantities that determine the fraction of the faint galaxies in rich clusters. We expect 
that decreasing 0:15 (i.e. steeper faint end slope) can increase the fraction of faint galaxies in 
rich clusters. It is therefore worth investigating if a more negative 0:15 than the value —1.32 
of the Ml model can save the halo model of Y03. 

It is important to emphasize that the conditional luminosity function method proposed 
by Yang et al. is a highly valuable approach to infer the galaxy distribution within dark 
halos from a wide range of observations. The discrepancy between the Ml model and our 
PVD luminosity dependence (especially the bimodal shape of the PVD) does not necessarily 
mean that the halo model approach is wrong or useless. Instead it does mean that the 
luminosity dependence of the PVD provides an important test for galaxy formation models 
that is independent of the correlation functions. Because the PVD is more sensitive to the 
small number of galaxies in rich clusters than the correlation function (Mo, Jing, & Borner 
1993; Mo, Jing, & Borner 1997), combining these two types of statistical quantities can put 
a more stringent constraint on the parametrization of the halo model. 

The galaxy-galaxy lensing of the SDSS constrains strongly the fraction of faint galaxies 
in rich clusters (Sheldon et al. 2004; McKay et al. 2002). The rapid decrease of the lensing 
signal with the luminosity of the lensing galaxies indicates that the average matter density 
around galaxies decreases with the decrease of galaxy luminosity (Tasitsiomi, et al. 2004; 
Seljak et al. 2004b). However, the galaxy-galaxy lensing observation may not be inconsistent 
with our PVD observation that quite a fraction of faint galaxies are in rich clusters. The 
main reason is that these two quantities as well as the two-point correlation function depend 
on the halo occupation number of galaxies differently. In other word, the faint population in 
rich clusters has a different weight on these different quantities. Because of its pair-square 
weighting nature, the PVD statistic is mostly sensitive to those faint galaxies in rich clusters 
(Mo, Jing, & Borner 1997). Because the mean lensing signal (Sheldon et al. 2004) is galaxy 
number weighted, the lensing statistic may be much less sensitive to the faint population 
in rich clusters than the PVD. Figure 6 also shows that the two-point correlation predicted 
by the current halo model (Ml model) is about 50 percent lower than the observed value 
for faint galaxies, indicating there is room for increasing faint galaxies in rich clusters. It 
would certainly be very interesting to examine quantitatively if such a halo model can be 
constructed to consistently interpret all these three quantities. Because these quantities 
depend on the halo occupation number in different manners, we expect that they will play 
important but complementary roles in constraining the halo occupation of galaxies. 
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6. Conclusions 

The analysis of the velocity fields of the galaxies in the 2dFGRS (COl) has led to 
a surprising discovery: The random velocities of the faint galaxies are very high, around 
700kms~^, reaching similar values as the bright galaxies. At intermediate luminosities the 
velocities exhibit a well defined steep minimum near 400kms~^. 

It seems that the galaxies in different luminosity intervals appear as different populations 
in their own right, defined by objective statistics. A look at Fig. (8) shows convincingly that 
this is actually the case. For this figure we have sorted the galaxies in 10 luminosity bins, 
each one magnitude wide, from magnitude —16.5 to —22 and plotted the value of (T12 at 
a wave number of A; ~ 1 /iMpc^^. Such a finely resolved binning of galaxies in samples of 
different luminosities is possible for the 2dFGRS redshift survey, because it is big enough to 
contain sufficiently many galaxies in each luminosity class. In each bin the PVD is a well 
defined quantity which can be measured reliably ( the very luminous galaxies have large 
error bars, because their are very few pairs of such objects at those scales). 

The PVD is an indicator of the depth of the local gravitational potential. Therefore 
we find the interesting result that the bright and the faint galaxies move in the strongest 
gravitational field. A substantial fraction of them must be in clusters, while most galaxies 
around magnitude M^ — blog^^h = —20.5, the M* galaxies in the Schechter luminosity 
function, populate the field. 

The bimodal nature of the correlation between PVD and luminosity may be used as a 
stringent test of galaxy formation models. We have investigated the halo occupation model 
(Y03 ) which has been optimally fitted to reproduce the luminosity function, and the two- 
point correlation function of the 2dFGRS. If we adapt this model to the PVD value of the M* 
galaxies, we see that it cannot give the high values found for the fainter galaxies. The PVD 
values of the model actually run opposite to the data and the model assigns smaller values 
to the fainter galaxies. This indicates that the assignment of galaxies to the dark matter 
halos must be done in a more intricate way as up to now. In Table 2 we have listed the 
PVD values for the different luminosity bins. These must be reproduced by an acceptable 
model for galaxy formation. The number of faint galaxies in clusters must be increased 
substantially to at least recover the high PVD found for them. One possible solution for 
the halo model is to raise the faint end slope of the conditional luminosity function in rich 
clusters. Also, the low value of 400 km s~^ found for the galaxies with magnitude —20.5 must 
mean that these galaxies reside in dark matter halos of galactic size. 

Another way, widely used, to connect dark matter to galaxies is the semi-analytic mod- 
eling (e.g., Kauffmann et al. 1999; Cole et al. 2000; Benson et al. 2000; Somerville & Primack 
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1999; Diaferio et al. 1999; Springel et al. 2001), where the dark matter distributions obtained 
from N-body simulations are supplemented with some of the physical processes important in 
galaxy formation using semi-analytic techniques. A test of the PVD vs luminosity for this 
type of models will be the aim of a subsequent paper. 

We have shown here in addition that our novel method of deriving the real space power 
spectrum and the PVD from the redshift space power spectrum allows to determine these 
quantities precisely and reliably from the 2dFGRS. The size of the 2dFGRS permits an 
investigation of the luminosity dependence of these quantities. New constraints on galaxy 
formation models can be derived from that. We have shown that the luminosity dependence 
of the clustering can be reproduced by an adequately chosen halo occupation model. This 
model fails, however, to reproduce the luminosity dependence of the PVD. 

We are grateful to Houjun Mo for a very stimulating and helpful discussion on the 
halo model, Uros Seljak for his helpful comment on the galaxy-galaxy lensing observation, 
and Xiaohu Yang for providing the halo model parameters of the low erg model and for 
some instructions of implementing the halo model in N-body simulations, and the anony- 
mous referee for helpful comments that improve our presentation. JYP would like to thank 
the Max-Planck Institute fiir Astrophysik for its warm hospitality during the time when 
this work was completed. GB wants to thank the Shanghai Astronomical Observatory for 
friendly infrastructure provided during his stay in Shanghai. The work is supported in 
part by NKBRSF (G19990754), by NSFC (Nos. 10125314, 10373012), and by the CAS-MPG 
exchange program. 
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Table 1: Samples selected according to luminosity 











South 


North 


Total 


Sample 


M,- 


- 5 logio h 




(no. of galaxies) 


(no. of galaxies) 


(no. of galaxies) 


1 


-16.5 > Mb - 


- 5 log^o h > 


-17.5 


5218 


3587 


8805 


2 


-17.0 > Mfc - 


- 5 logio /i > 


-18.0 


9314 


6078 


15392 


3 


-17.5 > Mb - 


- 5 logio h > 


-18.5 


14593 


11097 


25690 


4 


-18.0 > Mb - 


- 5 logio h > 


-19.0 


23703 


18670 


42373 


5 


-18.5 > Mb - 


- 5 logio h > 


-19.5 


34481 


25683 


60164 


6 


-19.0 > Mb - 


- 5 logio h > 


-20.0 


40995 


29241 


70236 


^7 

i 


-19.5 > Mk - 


- 5 logio h > 


-20.5 


40182 


27223 


67405 


8 


-20.0 > Mb - 


- 5 logio h > 


-21.0 


30934 


19204 


50138 


9 


-20.5 > Mb - 


- 5 logio ^ > 


-21.5 


15388 


8848 


24236 


10 


-21.0 > Mb - 


- 5 logio ^ > 


-22.0 


3739 


2162 


5901 



Table 2: The results of the PVD at A; = 1 /iMpc ^ of the luminosity subsamples 





Median Mag. 


a^[(3 = 0.45]^ 


a4/3(L)]b 


Sample 


Mb - 5 logio ^ 


(kms-i) 


(kms-i) 


1 


-17.15 


698 ± 103 


717 ±102 


2 


-17.62 


642 ± 73 


659 ± 72 


3 


-18.12 


723 ± 88 


737 ± 87 


4 


-18.61 


686 ± 42 


697 ± 42 


5 


-19.06 


609 ± 23 


616 ±23 


6 


-19.51 


498 ± 24 


500 ± 24 


7 


-19.96 


468 ± 21 


463 ± 21 


8 


-20.39 


424 d= 17 


412 ±17 


9 


-20.76 


599 ± 44 


580 ± 44 


10 


-21.18 


979 ± 303 


945 ± 296 



"The PVD determined with f3 = 0.45. 

''The PVD determined with the luminosity dependence of /3 taken into account. 
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Fig. 1. — Illustration of the smoothing effect on the determination of the redshift space power 
spectrum P{k, n) when the redshift space two-point correlation function ^(r^, tt) is weighted 
by a Gaussian function of width R = 20 /i"^Mpc. The solid lines, computed with Eq.(lO), are 
the P{k,fi) with the smoothing effect, compared to the dotted lines for P{k,fi) without the 
smoothing effect. In this plot, we use the linear CDM power spectrum of F = 0.2, (3 = 0.45 
and ay = 500kms~^ as input, k is in units of /iMpc^^ 
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Fig. 2. — The simulation test of the statistical method used in this paper that measures 
P{k,fi) from ^(rpjTr). The solid lines are the measurement of P{k,fi) for a simulation 
sample of galaxies with —19.5 < Mf, — Slog^Q/i < —18.5 based on the Gaussian weighted 
C{rp, tt) measurement, compared with the direct measurement of P{k, fi) based on the Fourier 
transformation of the galaxy density field in redshift space (the symbols). From top to 
bottom, the wavelength k is 0.10,0.16,0.25,0.40,0.63,1.0,1.6,2.5, and 5.0/iMpc"^ respec- 
tively. The galaxies are produced with the halo model, and the simulations have a boxsize 
L = 300 /i-^Mpc. 
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Fig. 3. — The luminosity function of galaxies generated with the halo model, compared with 
the observation of 2dFGRS. To the resolution limit Mj, — 5 log^g h = —18.5 for the simulation 
of L = 300 /i~^Mpc and M^ — 5 log^g ^ = —16.5 for the simulation of L = 100 h~^Mpc, the 
luminosity functions of the mock galaxies agree well with the observed one, and agree with 
each other in the two simulations. 
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Fig. 4. — The resolution effect on the predicted P{k,fi) in simulations of different boxsizes. 
Galaxies of —19.5 < M^ — 5 log^o h < —18.5 are analyzed. The results based on the 2dFGRS 
mock samples generated with L = 100 h~^Mpc simulations are plotted in symbols and those 
with L = 300 h~^Mpc simulations are plotted in connected lines. From top to bottom, the 
wavelength k is 0.10, 0.16, 0.25, 0.40, 0.63, 1.0, 1.6, 2.5, and 5.0 hMpc'^ respectively. 
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Fig. 5. — The redshift space power spectrum P{k, /i) measured in 2dFGRS. The symbols are 
for the whole survey, the dotted lines for the south subsample, and the dashed lines for the 
north subsample. The errors are plotted only for the whole survey that is estimated with 
the bootstrap method. The smooth solid lines are the best fits of Eq.(2) to data of the whole 
sample, (a) for -18.5 < Mb - blog^oh < -17.5; (b) for -19.5 < Mb - blog^oh < -18.5; 
(c) for -20.5 < Mb - 5 log^o /^ < -19.5; (d) for -21.5 < A'4 - blog^oh < -20.5. In each 
panel from top to bottom, the wavelength k is 0.2, 0.32, 0.50, 0.79, 1.26, 2.0, and 3.2 hMpc~^ 
respectively. 
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Fig. 6. — The real space power spectrum P{k) measured in the 2dFGRS (symbols): a) for 
-18.5 < Mfe - 51ogio/i < -17.5, b) for -19.5 < M^ - 5\ogh < -18.5, c) for -20.5 < 
Mf, — 51og^o ^ < —19.5, and d) for —21.5 < Mf, — 51og]^Q h < —20.5. The error bars of the 
observed results are given by the mock samples, as described in the text. The smooth lines 
are the predictions based on the halo model. 
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Fig. 7. — The PVD of galaxies measured in the 2dFGRS (symbols), compared with the 
predictions based on the halo model (solid lines). The error bars of the observed results are 
given by the mock samples, as described in the text. 
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Fig. 8. — The PVD measured at /c = 1 /iMpc~^ in the 2dFGRS (symbols for the whole sample, 
dotted line for the south, and dashed line for the north), compared with the predictions based 
on the halo model. The thick solid line is for the nominal model of ag = 0.9, the thick dashed 
one is for the model of o^ = 0.7, and the thick dotted one is for the model of ag = 0.7 but 
with Qq = 0.3. The error bars of the observed results are given by the mock samples, as 
described in the text. 
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Fig. 9. — The PVD measured at A; = 1 /iMpc ^ in the 2dFGRS. Symbols have the same 
meaning as in Fig. (8). Here the parameter /3 varies with luminosity as in Norberg et al. 
(2002a). 



